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ABSTRACT 



Aims. To search for extrasolar planets around intermediate-mass stars, we are conducting a precise RV survey around K-giants. 
Methods. We present high-accuracy RV measurements of a Ari made from November 2003 to February 2010. This star belongs to 
our sample of 55 K-giants studied for extrasolar planet and pulsation searches using the fiber-fed Bohyunsan Observatory Echelle 
Spectrograph (BOES) attached to the 1.8-m telescope at Bohyunsan Optical Astronomy Observatory (BOAO) in Korea. 
Results. We find a planetary companion with long-period and low-amplitude radial velocity (RV) variations in oscillating K2 III 
star a Ari (HD 12929). We do not find the correlation between RV variations and chromospheric activity indicators (Ca II H & K 
region, Ha line). The bisector analysis also shows that the bisector velocity span (BVS) does not show any obvious correlation with 
RV variation but has periodic component that may be attributed to the rotation of the star. If the RV variation is indeed caused by a 
planetary companion, an orbital solution with a period of P = 381 days, a semi-amplitude of K = 41 m s _1 , and an eccentricity of e = 
0.25 fits the data best. Assuming a possible stellar mass of M+ = 1.5 M Q , we estimate the minimum mass for the planetary companion 
of m 2 sin i = 1.8 M Jup with an orbital semi-major axis of 1.2 AU. 

Conclusions. Our finding of a likely exoplanet in a Ari supports searching for extrasolar planets around giant stars with multiperiodic 
oscillations. 
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1. Introduction 

Ever since the first extrasolar planet around a main-sequence star 
was discovered by Mayor & Queloz (1995), about 500 extrasolar 
planets have been found so far. Most of them are found around 
main-sequence stars using the precise RV method, but several 
planets have only been found around giant stars. 

Detection around intermediate-mass stars evolving toward 
the giant stage is easier than around early stage or main-sequence 
stars thanks to the availability of many sharp absorption lines for 
high-precision RV measurements. Thus makes giant stars suit- 
able targets for extrasolar planet detection with the RV method. 
Frink et al. (2002) discovered the first planetary companion 
around the K-giant star l Dra (K2 III), and thereafter, several 
companions around K-giant stars have been reported using the 
precise RV method (Setiawan et al. 2003a, 2003b; Mitchell et al. 
2003; Hatzes et al. 2005, 2006; Reffert et al. 2006; Johnson et al. 
2007, 2008; Dollinger et al. 2007, 2009; de Medrios et al. 2009; 
and Sato et al. 2007, 2008a, 2008b, 2010). However, in K-giants 
the velocity variations caused by planetary companions can be 
blended with the stellar pulsations and surface activities, which 
complicates identification of planetary companions. 

K-giant a Ari (HD 12929, HR 617) is one of our sample 55 
K-giant stars for which we have obtained precise RV measure- 
ments in the past eight years at BOAO in Korea. In previous pub- 
lication, we (Kim et al. 2006) reported multiperiodic oscillations 
in a Ari with periods P\ = 0.571 day (or aliases at 0.445 or 0.871 
days) observed during November to December 2003. The aver- 
age amplitude variation ranges ~ + 20 m s _1 just overnight. This 



variability can only be the result of stellar oscillations because 
of the short time scale. We continued to observe a Ari, by get- 
ting several spectra every night to search for long-term variation. 
Here we report the detection of a long-period and low-amplitude 
RV variation of a Ari possibly due to a planetary companion. 



2. Observations and analysis 

We acquired 678 spectra of a Ari from November 2003 to 
February 2010 using the fiber-fed high resolution (R = 90 000 at 
5000 A) echelle spectrograph BOES (Kim et al. 2007) attached 
to the 1.8-m telescope at BOAO in Korea. An iodine absorption 
cell (I2) was used to provide the precise RV measurements. Each 
estimated signal-to-noise ratio (S/N) at I2 wavelength region is 
about 250 with typical exposure time ranging between 90 and 
180 seconds. 

The extraction of normalized 1-D spectra was carried out 
using IRAF (Tody 1986) software. Further data processing, I2 
analysis and precise RV measurements were undertaken using a 
code called RVI2CELL (Han et al. 2007) developed at the Korea 
Astronomy & Space Science Institute (KASI). 

To demonstrate a long-term stability of the BOES we show 
relative RV measurements of the standard star r Ceti in Figure 
\T\ It shows that RV of r Ceti is constant with an rms scatter of 
6.7 m s _1 , over the time span of our observations. 
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Table 1. Stellar parameters for a Ari. 
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Fig. 1. RV measurements of a standard star r Ceti from 2003 to 
2010. 

3. The properties of a Ari 

We refer to the main photometric parameters of a Ari to Allende 
Prieto & Lambert (1999) taken from the Hipparcos catalogs of 
stars within 100 pc from the Sun and the main spectroscopic 
parameters to Hekker and Melendez (2007) taken at Lick obser- 
vatory. 

The mean RV of a Ari is - 14.6 km s _1 (Evans 1967). 
Berdyugina et al. (1991) derived r eff = 4850 + 150 K. And, Gray 
et al. (2002) determined r eff = 448 1 ± 1 1 K and [Fe/H] = - 0.221 
+ 0.036. Mozurkewich et al. (1991) measured an angular diam- 
eter for a Ari of 6.412 + 0.064 mas at 800 nm. The Hipparcos 
parallax of 49.48 + 0.99 mas results in a stellar radius of 13.9 
+ 0.3 R G . The most recent determinations of the projected rota- 
tional velocity v rot sin i = 3.44 km s _1 by Hekker & Melendez 
(2007) and 3.1 ± 1 km s" 1 by Carney et al. (2008) are close to 
each other. Based on the latter value that is given with error es- 
timation, we derived the range of upper limit for the rotational 
period of 

P ro , = 27rR / (v ro t sm = 172 - 335 days. 

Berdyugina et al. (1991) determined log g = 2.6 + 0.3, almost in 
agreement with log g = 2.57 + 0.2 as determined by McWilliam 
(1990). Surface gravity and its uncertainty log g = 2.6 + 0.3 
yields a stellar mass in the range of 1.4 - 5.6 M G . To estimate 
the mass of a Ari independently, we determined atmospheric 
parameters directly from our spectra. By using measured 268 
equivalent widths (EW) of Fe I and Fe II lines, we estimated 
r e ff, [Fe/H] = - 0.09 + 0.03, log g, and v, mc , of the star using 
the program TGVIT (Takeda et al. 2005). Our estimated surface 
gravity with its formal uncertainty is log g = 2.33 + 0.06. That 
yields a stellar mass of 1.5 + 0.2 M . Table 1 summarizes the 
basic stellar parameters of a Ari. 

4. Orbital solution 

The long-term variation of a Ari was first noted by Walker et 
al. (1989) and Larson et al. (1999) who used CFHT 3.6-m and 
DAO telescopes with a hydrogen-fluoride (HF) absorption cell. 
These data were insufficient (45 spectra) to draw any conclusions 
about possible periodic variability. We recalculated the standard 
deviation of their 45 measurements of RV variation of o~ cfht = 
41.6 + 6.2 m s during 13 years. The standard deviation of our 
678 measurements spanning eight years is <t BO ao - 27.1 ±1.6 



Parameter 


Value 


Reference 


Spectral type 


K2 III 


Hekker & Melendez (2007) 


m» [mag] 


2.01 


Hipparcos 


M v [mag] 


0.48 ± 0.09 


Hipparcos 


5-V[mag] 


1.151 


Hipparcos 


Parallax [mas] 


49.48 ± 0.99 


Hipparcos 


Diameter [mas] 


6.412 ± 0.064 


Hipparcos 


Distance [pc] 


20.2 


Hipparcos 


RV [km s" 1 ] 


-14.6 ±0.1 


Evans (1967) 


r eff [K] 


4480 


McWilliam (1990) 




4850 ± 150 


Berdyugina et al. (1991) 




4481 ± 11 


Gray et al. (2002) 




4553 ± 15 


This work 


[Fe/H] 


-0.25 


McWilliam (1990) 




-0.22 


Mozurkewich et al. (1991) 




- 0.09 ± 0.03 


This work 




2.57 ±0.2 


McWilliam (1990) 




2.6 ± 0.3 


Berdyugina et al. (1991) 




2.33 ± 0.06 


This work 


Af* [Ma] 


1.5 ±0.2 


This work 


R+ [Re] 


13.9 ±0.3 


derived 


v rot sin i [km s -1 ] 


3.44 


Hekker & Melendez (2007) 




3.1 ± 1 


Carney et al. (2008) 


P rot 1 sin i [days] 


172-335 


derived 


v micm [km s -1 ] 


1.7 


McWilliam (1990) 




1.40 ± 0.06 


This work 



Figure[3]shows the Lomb-Scargle periodogram (Lomb 1976; 
Scargle 1982) of the original RV time series for a Ari. It shows 
a significant power at f { =0.00263 cd" 1 (P = 380.8 days). The 
individual RV measurements are shown on the top panel in 
Figure|2] and they exhibit clear variations. 

The statistical significance of this signal was calculated 
using the bootstrap randomization technique (e.g. Hatzes & 
Mkrtichian 2004 and references herein). The measured RV val- 
ues were randomly shuffled keeping the observed times fixed. 
The Lomb-Scargle periodogram was then computed 200 000 
times for each "random" data set. The fraction of a periodogram 
having power higher than the amplitude peak in the range of < 
/ < 0.03 c d -1 represents the false-alarm probability (FAP) that 
noise would create the detected signal. A small FAP is therefore 
an indication of a significant peak. We found that among 200 000 
trials the highest peak of periodograms never exceed a signal in 
the original data, thus the FAP < 5 x 10~ 6 . 

The solid line in the top panel of Figure [2] is the best-fit 
Keplerian orbit with P = 380.8 ± 0.3 days, semi-amplitude K = 
41.1 + 0.8 m s _1 and eccentricity of 0.25 ± 0.03. We estimated 
the minimum mass of the companion to be 1.8 + 0.2 Mj up at a 
distance of 1 .2 AU from a Ari. The RV phase diagram for the or- 
bit in Figure|2]is shown in Figure|4] CFHT data (open circles) are 
also shown in the phase diagram to the period of BOES data set. 
Uncertainties of CFHT RV measurements are 2 ~ 3 times larger 
than that of BOES and also do not conform perfectly with BOES 
phase diagram because of low instrumental exquisiteness and 
some noise pattern of the data (Larson et al. 1993, 1999). The 
rms residual from the orbital solution is 17.8 m s _I (Figure 0. 
This is almost certainly due to the short-period variations and 
intrinsic variability of stellar oscillations, as shown by Kim et al. 
(2006). All the orbital elements are listed in Table 2. 
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Fig. 2. RV measurements (Top) and rms scatter of the residuals (Bottom) for a Ari taken at BOAO from November 2003 to February 
2010. The solid line is the orbital solution with a period of 380.8 days and an eccentricity of 0.25. 
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Fig. 3. The Lomb-Scargle periodogram of the RV measurements 
for a Ari. The periogram shows a significant power at a fre- 
quency of 0.00263 c d _I corresponding to a period of 380.8 days. 
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Fig. 4. RVs of a Ari phased to the period of 380.8 days. Filled 
and open circles denote BOES and CFHT spectroscopic mea- 
surements, respectively. The solid line is the orbital solution that 
fits the data with an rms of 17.8 m s _1 . 



5. The nature of the RV variations 

K-giant stars exhibit pulsation, as well as surface activity, re- 
sulting in low-amplitude RV variability on different time scales. 
While short-term period (hours ~ days) RV variations have been 
known to be the result of stellar pulsations (Hatzes & Cochran 
1998), long-term period (a hundred of days) RV variations with 
a low-amplitude may be caused by three kinds of phenomena: 
planetary companion, rotational modulation of surface features, 
or non-radial oscillation modes. To establish whether a low am- 
plitude and long-period RV variations are caused by surface ac- 
tivities or companions we examined: 1) the Hipparcos photom- 
etry, 2) the stellar chromospheric activity indicators like the EW 
of sensitive lines, and 3) spectral line bisectors measurement. 



Table 2. Orbital parameters of the RV measurements for a Ari 
b. 



Parameter 


Value 


Period [days] 


380.8 ± 0.3 


^periastron [JD] 


2451213.52 ± 7.36 


K [m s" 1 ] 


41.1 ±0.8 


e 


0.25 ± 0.03 


oj [deg] 


89.16 ±6.84 


f(m) [Mq] 


(2.476) x 10-" 


m sin i [M Jup ] 


1.8 ±0.2 


a [AU] 


1.2 


o-(O-C) [m s-»] 


17.8 



5.1. Hipparcos photometry 

We analyzed the Hipparcos photometry for a Ari in order to 
search for possible brightness variations due to the rotational 
modulation of cool surface spots in case the RV variation is 
caused by the rotation. For 2.5 years, the Hipparcos satellite ob- 
tained 66 photometric measurements of a Ari, and we used the 
whole data except one measurement with a large error. 



Figure[5]shows the Lomb-Scargle periodogram of Hipparcos 
photometric measurements calculated in the frequency interval 
of - 0.03 c d -1 . There are no significant peaks near the fre- 
quency of f\ = 0.00263 c d -1 that corresponds to the period 
of 381 days. Peaks in the periodogram are the typical low- 
frequency noise peaks so they cannot be considered as signifi- 
cant. 
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Fig. 5. The Lomb-Scargle periodogram of the Hipparcos photo- 
metric measurements of a Ari for 2.5 years. The vertical dotted 
line marks the location of the orbital frequency. 



5.2. Chromosphehc activity indicators 

The EW variations of Ca II H & K, H a , and Ca II 6882 A 
lines are frequently used for chromospheric activity analysis. 
The emissions in the Ca II H & K core are formed in the chro- 
mosphere and show a typical central reversal in the existence of 
chromospheric activity (Pasquini et al. 1988; Saar & Donahue 
1997). Because the Ca II H & K line region of the BOES spectra 
does not have enough S/Ns to estimate EW variations of each 
spectrum reliably, we used the averaged Ca II H & K line pro- 
files at the upper and lower part of RV curve to see any system- 
atic difference related to the RV variations. We selected spec- 
tra for a Ari that correspond to the positive deviation on JD- 
2453779.9451 (RV = 52.9 m s -1 ) and the negative deviation on 
JD-2454027.1307 (RV = -48.1 m s" 1 ) of RV curve. These spec- 
tra in the Ca II H & K region were plotted in the Figure [6] Core 
reversal in the Ca II H & K central region is not obvious and 
there is no systematic difference in the line profiles at these two 
extremal parts of the RV curve. It means that K-giant a Ari, at 
the moment of observations, exhibits a modest chromospheric 
activity. 

The Balmer H a line can also be a useful indicator of chromo- 
spheric activity to determine whether the RV variation is due to 
intrinsic stellar variability (Kiirster et al. 2003). Unlike in main 
sequence stars where line is blended with telluric and broad 
lines, it is easy to measure the EW of H a line in giant stars. The 
mean EW of the H a line in a Ari is measured to be 1087.5 ±9.1 
mA. The rms of 9. 1 mA corresponds to a 0.84 % variation in the 
EW. Figure [7] shows the EW of the H„ line as a function of RV. 
A gradual slope (solid line) of the fitted line shows that the EWs 
of the H a line are not correlated with RV variations. 

We conclude that the chromospheric activity and surface 
spots, if present, have too small amplitude to be resolved in the 
Call H & K and the H„ line EW variations. Another activity in- 
dicator Ca II 8662 A line is not suitable because the significant 
fringing and saturation of our CCD spectra at wavelengths are 
longer than 7500 A. 

5.3. Line bisector variations 

Queloz et al. (2001) report that the analysis of the line bisector 
can distinguish the origin of the observed RV variations, which 
is very useful for the messy K-giant stars. The RV variations 
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Fig. 6. The spectral region of the Ca II H & K lines in a Ari. 
Two spectra exhibiting the positive (top) and the negative (mid- 
dle) RV deviations are displayed (see text). Bottom shows the 
differential spectrum. Any systematic difference in line profiles 
and core reversals are not apparent around the Ca II H & K cen- 
tral region, indicating modest activity. 
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Fig. 7. Relative EW(A) vs. RV of the H a line for a Ari from 
September 2004 to February 2010. The diagram shows no cor- 
relation between the EWs and the RVs. 



due to planetary companions should not produce any changes in 
the spectral line shape, whereas the surface inhomogeneities do. 
Stellar rotation modulations of surface features can create vari- 
able asymmetries in the spectral line profiles. Thus, the analysis 
of the shapes of spectral lines may prove or disprove the exis- 
tence of planetary companions. 

Bisector velocity span (BVS) is defined as the difference 
between two regions of the bisector at the top and bottom of 
the profile to measure changes with high accuracy. To search 
for variations in the spectral line shapes, we measured the BVS 
using the mean profile of the spectral lines. We calculated the 
mean profile of spectral lines using the least squares deconvo- 
lution (LSD) technique (Glazunova et al. 2008). We used the 
Vienna Atomic Line Database (VALD; Piskunov et al. 1995) to 
prepare the list of spectral lines. The total of 4339 lines within 
the wavelength region of 4500 - 4900 and 6450 - 6840 A were 
used to construct the LSD profile, which excluded spectral re- 
gions around the I2 absorption region, hydrogen lines, and re- 
gions with strong contamination by terrestrial atmospheric lines. 
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Fig. 8. BVS variations for a Ari from November 2003 to 
February 2010. The solid line shows a sine wave with a period 
of 211 days fit. 



We estimated the BVS of the mean profile between two different 
flux levels of central depth levels 0.8 and 0.25 as the span points. 

The BVS vs JD shown in Figure [8] exhibits some periodic 
variations. The BVSs of the JD 2454453 - 2454506 spectra seem 
to be shifted with respect to the preceding and subsequent data 
by the amount grater than the + 100 m s -1 . We suspect that this 
bias was come from the temporary mis-alignments of the BOES 
echelle spectrometer. But, we did not find any BVS shift in the 
K-giant -y 1 Leo for the same observing epoch as where a plane- 
tary companion system was detected by using the same echelle 
spectrometer BOES (Han et al. 2010). The RV code RVI2CELL 
takes any changes into account in the instrumental profile in 
calculations of RVs, and thus instrumental misalignment is not 
likely to have affected the RV data. 

To search for periodic signals in the BVS data we calcu- 
lated the Lomb-Scargle periodogram of the BVS (top panel in 
Figure [5). The largest peak in the Lomb-Scargle periodogram is 
at a very low frequency, fa = 0.00108 c d" 1 (P = 926.8 days) 
with an FAP less than 9 x 10~ 6 . After removal of this signal, 
we plot the Lomb-Scargle periodogram of the residuals at bot- 
tom panel in Figure [9] The highest peak at fa - 0.00474 c d _I 
(P = 211.2 days) has high FAP ~ 2 x 10~ 5 . Figure |8] shows the 
individual fit of the frequency to BVS data. 

Taking the stellar radius and rotational velocity (y rol sin i) 
value into account we derived the expected range of rotational 
period of 172 - 335 days. Found about 211 day periods of BVS 
variations well fit this interval. We can summarize that the vari- 
ation can be attributed to stellar surface active regions, inhomo- 
geneities, and rotation. 

6. Discussion and conclusion 

Our analysis of the RV measurements for a Ari revealed a 381 
day periodic variation that persisted between 2003 - 2010 for 
almost six cycles with no changes in phase or amplitude. We 
also confirmed the short-period oscillations previously found in 
a Ari (Kim et al. 2006). We found no evidence of the correla- 
tion between the 381 day RV variation and the EWs of H a lines 
or the Hipparcos photometry. But, we found about 211 day pe- 
riodic variations in the BVS that we attribute to the rotational 
modulation of line profiles from the surface spots and about 927 
day variations in the unknown origin. Thus, we conclude that the 
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Fig. 9. The Lomb-Scargle periodogram of the BVS variations 
from November 2003 to February 2010 (top panel) and after 
subtracting the low-frequency (fa = 0.00108 c d _1 ) variations 
(bottom panel). The highest frequency is at fa = 0.00474 c d -1 
(P = 21 1.2 days). The vertical dotted line marks the location of 
the orbital frequency. 

planetary companion in a Ari is the most likely cause of the 381 
day RV variations. 

If confirmed, the determined orbit of the 1.8 Mj up mass 
planet is fairly eccentric, e = 0.25 ± 0.03, but, up to now, 6 out of 
22 companions around K-giants show more than 0.2 eccentric- 
ity: e = 0.7 of i Dra (Frink et al. 2002), e = 0.2 of Kl III star HD 
47536 (Setiawan et al. 2003), e = 0.27 of Kl II-III star HD 13189 
(Hatzes et al. 2005), e = 0.43 of Kl III star 4 UMa (Dollinger et 
al. 2007), e = 0.38 of K1.5 III star 42 Dra (Dollinger et al. 2009), 
and e = 0.46 of K2 III star HD 1 10014 (de Medeiros et al. 2009). 
The eccentricity for the planetary orbit of a Ari is therefore not 
unusual among K-giants, and such eccentric modulation is diffi- 
cult to produce by surface variations or stellar oscillations. 

The source a Ari belongs to yet smaller group of planet 
hosting K-giants showing short-period oscillations (e.g. ySGem 
(Hatzes & Zechmeister 2007), l Dra (Zechmeister et al. 2008)). 
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